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In oflfering to the student a treatment of the Problem of 
the Earth's rotation somewhat diflferent from that which 
has been usually given in elementary text-books, «, few 
words of explanation . are necessary. 

The FIRST PART consists of an application of the method 
of the variation of elements to the general problem of 
rotation. That the formulae for calculating these elements 
are identical in the motions both of translation and rotation, 
appeared so remarkable, that it might be well to present 
the latter in a form easily accessible. As far as I am aware, 
an elementary investigation of these formulae has not yet 
been given : in attempting to supply this, I have adopted 
a method somewhat similar to that which I have given for 
the corresponding equations of the motion of translation, in 
an Elevientary Treatise on the Planetary Theory. The striking 
analogy, thus developed, between the solutions of problems, 
in appearance so dissimilar, may, I hope, lead to a more 
complete study of Lagrange's beautiful theory of the varia- 
tion of arbitrary constants. 

In the SECOND part the general rotation-formulae are 
applied to the particular case of the Earth. These formulae 



VI PREFACE. 

afford a simple and accurate proof of the important theorem 
of the Stability of the Earth's axis and of the motion about 
it^ so far as these depend upon the attractions of distant 
bodies. In this I have followed M. de Pont^coulant, The 
remaining pages are devoted to a consideration of the motion 
of the Earth's axis in space. In this I have obtained the 
formulae for calculating Precession and Nutation, first, by 
an appUcation of the general method, and afterwaa-ds. by 
an independent process ; but I have not carried the approxi- 
mation further than has been usual in elementary text- 
books. 

C. H. H. CHEYNE. 



I, Dean's Yabd, 
Wbstminbtbb, 
September, 1867. 
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PART L 

aENERAL THEORY OP ROTATION* 

1. In that part of Physical Astronomy which usually 
goes by the name of the Planetary Theory we are con- 
cerned with the motions of translation only of the planets 
in space : we now propose to consider their motions of rota- 
tion. The principles of the conservation of the motions of 
translation and rotation permit us to consider these isepa- 
rately, and to treat the latter as if the centre of gravity were 
a fixed point. We shall adopt a method perfectly rigorous, 
and free from all assumptions, with the single exception 
of the hypothesis, already required in the Planetary Theory, 
that the. attracting bodies are so distant that their action 
may be supposed the same as it would be if their whole 
masses were condensed into their centres of gravity. Thus 
we shall obtain, for the determination of the motion, for- 
mulae applicable to the case of any planet or other rigid 
body : an interesting application of these will then be afforded 
by the special circumstances which occur in the EartKs 
Motion of BotcUion. 

2. If the planets were exactly spherical in shape, it is 
clear that the attractions of the Sun, Moon, and of the other 
planets could produce no effect upon their rotation, since 
they would all pass through the centre of gravity. But 

c. ^ 
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although this is not the case, yet the deviation from exact 
sphericity being very small, the motion will dififer only 
slightly from what it would be if these disturbing forces did 
not exist. We shall, therefore, by neglecting. them obtain 
first an approximate solution of the problem, and then by 
the method of the variation of parameters deduce from it 
accurate results. Since, however, the motion of a rigid body 
about a fixed point under the action of no forces is dis- 
cussed in works on Rigid Dynamics, we shall here consider it 
only so far as is necessary for the purpose of obtaining re- 
sults which will be required in the sequel. 



Undisturbed Motion^ 

3. Let »i, «j, »3 be the angular velocities of a planet 
about the principal axes at its centre of gravity ; A, B, C 
the moments of inertia about these axes : then Euler's 
equations give 

5^-(C-^)a,.«, = 0, 

Multiplying these equations by a)^, co^, ©3 respectively, 
adding, and integrating, we have 

where h is the constant of integration. 
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Again, multiplying by Aod^, Bao^y God^, adding, and iute- 
grating, we have 

where i? is the constant of integration. 

From these two equations we obtain 

, T^ - Bh + {B - C) C(o* 
"» ' (,A-B)A 

t_ -lt^ + Ah + iC~A)Co>* 
**« ~ {A-£)B 

Substituting these values, the third equation of motion 
becomes 

^<&»,_ {y-m + (B-C) C<o,*}*\-k^+Ah + {C-A)Ca,'}* 
dt~ ^{AB) 

whence t + l 

where I is the constant of integration. 

This integral cannot in the general case be found; we 
may however approximate : thus t is known in terms of a,, 
and consequently cd, in terms of t ; and then from above, 
cDj, a>. are also known. 

4. With respect to the constants introduced by the inte- 
gration, we may remark that h represents the vis viva 
(Routh's Bigid Dynamics, Art. 194), and k the area conserved 
on the invariable plane. To prove the latter point, the 
areas conserved on the principal planes being Aa)^, B(o^, Ca>^ 
(Routh's jS. J)., Art 179), and the direction cosines of the 
invariable plane with reference to the principal axes 

-4q)j B(o^ Cto^ 

i6 IC iC 
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(Routh's iJ. D., Art, 125), the area conserved on the invari- 
able plane 

= ^ -| r— ^ -f- " .. 

5, When cd^, o),, o), are known at any time, the result- 
ant angular velocity of the planet is known, and also the 
position of the instantaneous axis of rotation with reference 
to the principal axes. It remains to shew how the position 
pf these axes in space may be determined. 

Suppose a sphere described with its centre at the centre 

of gravity of the planet and its radius of any magnitude : 

take as a plane of reference any fixed plane passing through 

the centre of gravity, and let it cut the sphere in the great 

circle ON; also let the principal plane of xy cut the sphere 

in the great circle NAB, -^T being the node of this plane 

upon the fixed plane, and A, B, C the points where the 

sphere is cut by the principal axes of x, y, z. Take P the 

pole of ON, and join FA, PB, CA, OB by arcs of great 

circles, 

P 

O 




Let the angle ONA^ 0, 0N= f, NA = <I>: then if the 
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angles By ^, y^ be known, the position of the planet will be 
determined with reference to the fixed plane. 

Now we may consider the planet to be moving with angu- 
lar yelocities. a>^, o),, cd, about the principal axes ; or with 

angular velocities -^ > 777 ' ;^ ' ^^^ ^^* about a normal to 

the fixed plane, the second about the line of nodes, the third 
about the principal axis of z. We shall adopt the usual 
convention with respect to signs, and consider positive those 
angular velocities which tend to turn the planet round the 
axes of X, y, z from y io z^z to a?, x to y, respectively. 

Thus, resolving about the principal axes, we have 



do • , rf'Jr T1T» 

0), s= -^ sm ^ — -^ COS PBy 



dt ^ dt 



d<b d'dt ^ 



Now, by Spherical Trigonometry, 

COS PA = — sin ^ sin <^, 
cos PB = — sin ^ cos ^ ; 

therefore a)j = — -^ cos ^ + -t^ sin tf sin ^, 

dO , t , dylr . ^ , 
», = -1- sm 9 + -^ sm 6 cos 9, 

d<b d^it A 
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Hence also, 



-^sin^ = a)jSin<^ + ©,cos <f), 

d0 

— s — 0)^ COS <^ + fi>, sin <f>^ 

dff> d^ ^ 



• ••*«• \^) • 



By substituting the values of co^ m^, ©3 obtained as above 
(Art 3), and then integrating these equations, 6, 0, and '^ 
would be determined, and thus the position of the principal 
axes at any time would be known. The integration, how- 
ever, canno in general be eflFected; so that we are obliged 
to have recourse to a special hypothesis with regard to the 
position of the fixed plane of reference. If we take for this 
purpose the invariable plane, the process becomes much sim- 
plified. 

6. Let then 0^, ^j, -^^ denote relatively to the inva- 
riable plane the same angles which relatively to the original 
plane of reference have been denoted by 6, ^, '^. Then, the 
direction cosines of the invariable plane with reference to the 
principal axes being 

Aw^ Boo^ do, 

rC ■ tC ' K 

respectively, we have (see figure of preceding Article) 

-jT^ = cos PA = — sin 0^ sin <f>^, 

-=-^ = cos PJ?= — sin ^jCOS <f>^^ 
-^ = cos^,; 

• There is much diaagreement between writen as to the measurement of 
the angles employed in these kinematical equations ; the above, however, agrees 
with La Place, Poisson, and Pont^ooulant 
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therefore 

cos 0, =s -r-» , 

These equations give 0, and ^, : to obtain -^^ substitute 
in the first of the equations (^) of Art. 5 ; thus 

-^ sin* 0, = a>, sin 0, sin ^^ + ©, sin 0, cos ^,; 

therefore ^ (¥ - C*m^ = - {Aw* + 5«,*) * 

{h-C<o,*)k', 

combining this with the result of Art. 3, and integrating, we 
have 

where g is the constant of integration. 

Since cd, is known in terms of t from Art. 3, these equa* 
tions give O^y <^j, '^^\ so that the position of the principal 
axes is known at any time with reference to the invariable 
plane. Since, however, when the disturbing forces are 
taken into account, this plane ceases to be absolutely invari- 
able, it will be convenient to be able to refer the motion to 
some other plane which does remain fixed, and which may 
be taken as a plane of reference : this we can now do by 
Spherical Trigonometry. 



8* 
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7. Let the surface of a sphere of any radius, with its 
centre at the centre of gravity of the planet, be cut by the 
fixed plane of reference, the invariable plane, and the princi- 
pal plane of xy, in the great circles OMN, MI, IN A respec- 
tively. 




As before let ON^ '^, NA = <^, the angle ON A = ; also 
take M as the origin from which '^^ is measured, and let 
MI=^'^,, lA =0j, the angle MIN^ 0^ : let OM (the longi- 
tude of the node of the invariable plane) = a, and IMN (its 
inclination to the plane of reference) = 7. 

Then the sides of the spherical triangle IMN are 

and the angles respectively opposite to these, 

0xf w--^, 7. 

Hence, by the formulae of Spherical Trigonometry, we 
have 

cos — COS7 cos 0^ — sin 7 sin 0^ cos -^j, 

sin (<^j — <^) sin ^ = sin 7 sin '^1 , 

sin [^ — a) sin = sin 0^ sin -^i, 

which determine 0, ^, -^^ when 01^ j>^, -^j, a, 7 are known. 
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Disturbed Motion, 

8. Having now shewn how to determine the position 
and velocity of rotation of the planet on the hypothesis that 
no forces act upon it, we proceed to a rigorous treatment of 
the problem. We shall employ the principle of the varia- 
tion of parameters, and suppose the results already obtained 
to represent the true solution, the arbitrary constants or 
elements being no longer constants, but variable quantities, 
which it will be our object to determine. We shall speak of 
the forces which produce this variation as disturbing forces. 

The elements which have been already employed in the 
undisturbed motion are six in number, viz. h, k, 1, g, a,y: 
in considering these as variable we shall arrive at the very 
remarkable result that the equations for calculating their 
variations are precisely the same as the corresponding equa- 
tions for the motion of translation in the Planetary Theory, 

Def. The plane of which the direction cosines are 

-4a)j Bw^ C(o^ 

fC iC rC 

the area conserved upon which has been shewn (Art. 4) to be 
equal to k, will in future be termed the pkme of maximum 
areas, on account of the property which it possesses, that k 
is a maximum*^ since it can no longer be considered in- 
yariable. 

* See Bouth's Itigid Dynamia^ Art. 174* 
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9. To .find an expression for the sum of the moments of 
ihe disturbing forces about any line through the centre of 
gravity of the planet. 

We shall suppose the disturbing body (which may be the 
Sun, Moon, or another planet), so distant that it may be 
considered to attract as if condensed into its centre of gravity. 

Let m be the mass of the disturbed, m' of the disturbing 
body, pj the distance of the centre of gravity of the latter 

from an element Sm. of the former : also let F- = — . Then, 

Pi 
if a- denote the length of the arc of any curve measured 

from some fixed point to the element Bm^, the disturbing 

ferce on this element in direction of the tangent to the curve, 

and tending to increase a, will be 

If we suppose this arbitrary curve so drawn that its 
tangent at the point where the element is situated is perpen- 
dicular to the axis about which the moments are to be taken, 
and denote by p the distance of Sm^ from this axis, the 
moment of the force will be 

Sm,.p -^ . 

Let the small arc Scr subtend an angle hx ^t the nearest 
point of the axis ; then Scr =^S;^ and the moment becomes 

Similarly, if F, refer to an element Sm,, the moment of 
the disturbing force on this element will be 



Sm 






GENERAL THEORY OF ROTATION. 11 

Sx being the same as for the element Bm^ since these ele- 
ments are supposed rigidly connected. 

Hence the suto of the moments of the disturbing forces 
on all the elements of the planet 

or, if we write V for 2 (Sm,.Fj), the sum of the moments 
will be —J— , 

where F' = m'% — . 

P 

Cor. If there are several disturbing bodies m', m", &c. 
and V, V", &c. are the functions corresponding, the sum of 
the moments of the forces due to their action will still be 

^, where V=r+V" + ... 
dx 

The result of this Article may be thus enunciated : — Sup- 
pose a small arbitrary rotation given to the planet about any 
aats through an am,gle B^; then, supposing V expressed in 
terms ofx <^^ quantities which do not vary in this hypotheti- 
cal motion, the sum of the moments of the disturbing foroes 
abovit this axis will be expressed by the partial differential 

coej/icient -^ — . 

The function Fis thus clearly analogous to the disturbing 
function R of the Planetary Theory. 

10. "We may express F in various ways which will be 
foimd useful : 

(i) As a function of 0, up, y^. Let x, y, z he the co- 
ordinates of an element Sm of the disturbed planet, the fixed 



12 THE earth's motion OF ROTATION. 

plane of reference being that of xy^ and the axis of x thd 
line from which '^ and a are measured ; let x, y\ z be the 
co-ordinates of the centre of gravity of the disturbing body 
referred to the same axes; and let a?,, y^, «, be the co-ordi- 
nates of hm referred to the principal axes of the planet. 
Then 

F=w'2— , (Art. 9), 

r 



Now if \, fi, V be the angles which the fixed axis of x 
makes with the principal axes of a?, y, z, we have 

a? = ajj cos \ + y^ cos /A H- «j cos I' ; 

and by Spherical Trigonom,etry (see fig. to Art. 5), 

cos \ =s cos <^ cos -^ + sin ^ sin -^ cos 0, 
cos /A = — sin ^ cos -^ + cos ^ sin '^ cos 0, 
cos V = sin '^ sin 0^ 

Thus X may be expressed in terms of 0, <f>, y^ and of thc> 
co-ordinates of Bm referred to the principal axes. Similarly, 
y and z may be expressed in terms of 0, <f>, y^» If the values 
of X, y, z so obtained be now substituted in the expression 
for V, it will become a function of 0, <f>,'^; and, in so far as 
it depends upon the disturbed planet, of quantities indepen- 
dent of the time. 

(ii) As a function of 0^, <f>^, '^^, a, 7. By Art. 7, we 
have,. ~ ^ 

'f-a=/(^,, f„7), 
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the symbol/ denotiDg a diflFerent function in each case. If 
then we suppose V to have been expressed as a function of 
0, <f>, '^ by (i), these equations will enable us to express it as 
a function of 0^, <f>^, -^j, a, y. 

(iii) As a function of t and the elements. Collecting 
together the results of Arts. 6 and 3, and making Oo), =^, 
we may write 

^, =/(*,*), 

ir^ + ff='f{h,k,8), 

t+l=:f{h,kj8). 

Supposing Fto have been expressed as a function of 

^i> ^i» '^i> *' 7' ^y (^^)» ^^^ ^^^ three of these equations will 
enable us to express it as a function of a, h, k, g, a,y: then, 
if « be eliminated by means of the fourth equation, it will 
become a function of ^ + ?, A, k, g, a,yi that is, of t and the 
elements. 

11. We now proceed to obtain equations for calculating 
the values of the elements at any time, commencing with h, 
the element of vis viva. 

Let T denote the vis viva due to the rotation of the 

dV 
planet, Bm^ -— the resolved part of the disturbing force on 

an element Sm^ of its mass in the direction of motion of the 
element : then the equation of vis viva gives 

where in \ ^ , the differential coefficient is taken only in 
at 
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80 fax as Hs involved through the co-ordinates of the element 
Swij of the disturbed planet. This equation may be written 

if, as in Art. 9, we write V^^iJim^.V^. Now from the 
result of (iii) in the preceding Article we notice that t-^l 
always occurs in F as one quantity ; therefore 

d{V) _dV 
dt '^ dl'' 

also in the undisturbed motion r=A; hence our equation 
becomes 

d_h_ dV 

dt" dl' 

from which h may be calculated. 

12. This equation may also be written 

dh^^djV) 
dt dt ' 

under which form it is easily seen to be identical with the 
corresponding equation 



dt\ a) dt ' 



of the Planetary Theory. (See Planetary Theory^ Art. 26.) 
In fact it appears that the element of vis viva will be given 
by a similar equation in all cases of motion, whether of 
translation or rotation, when the disturbing bodies attract 
according to any law expressed by a function of the distance. 
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Eqyuxtiona of Motion, 

13. Let Aj, Ag, A3 denote the areas conserved about three 
rectangular axes of x, y, «, originating in the centre of 
gravity of the disturbed planet, and moving with angular 
velocities )8j, )8,, ^S, about their instantaneous positions ; also 
let L, My N be the moments of the disturbing forces about 
these axes: then (Eouth's Rigid Dynamics, Art. 120) we 
have the equations of motion 



dt 

dk, 
dt 



s 



-AA+^A=-^- 



Now let A3 denote the area conserved on the plane of 

maximum areas : then from the undisturbed motion (Art. 4) 

we have 

Ajg = A; 5 

and since the areas conserved upon all planes perpendicular 
to the plane of maximum areas are zero*, we must also have 

A, = 0, A, = 0, 

always ; and therefore 

* This follows from the fact that the area conserved upon any plane is equal 
to the projection upon it of the area conserved upon tfie plane of maximum 
areas. (Routh's Bigid Dyrmrnks^ Art. 174.) 
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hence our equations become 

dt • 

J 7 

• We. may express fi^ and fi^ in terms of -j- and ^ , a and 

7 being respectively the longitude of the node and the 
inclination of the plane of maximum areas to the fixed 
plane of reference. Since no assumption has yet been made 
Tvith regard to the position of the axis of x, it will be con- 
venient to suppose it to coincide with the line of nodes. 
Employing the equations of Art. 5, and making w^=s/3 , 
6>a = iS,» ^ = 7> ^ = 0, '^ = a, we have 



da . 
^sm7 = ^., 

dy 
It 



^ = -/3. 



and the equations of motion become 

A?-7r sm7 = 2y, 
at 

14. With respect to the last of these equations, we may 
remark that it is of the same form as it would be if the 
plane of maximum areas were still invariable. Now on 
referring to Art. 18 of the Planetary Theory it will be seen 
that the equations for the motion of translation, when re- 
ferred to axes in the plane of the orbit, take the same forms 
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as if this plane were at rest, and moreover, that the second 
equation there obtained is identical with the third of thfe 
preceding Article*, This coincidence is not accidental, but 
arises from the fact that in the motion of translatioQ, the 
plane of maximum areas is the plane of the orbit. In fact^ 
the equations obtained in the preceding Article for the motion 
of rotation are equally applicable to the motion of transla- 
tion, and we shall see that the resulting formulsB for the cal^ 
culation of the elements involved are identical in the two 
motions, 

15. It now only remains to find expressions for i, M, N 
in terms of the diflFerential coeflScients of the function V, in 
order to obtain from the equations of Art. 13 formulae for 
calculating the elements a, 7 and k. We shall suppose V 
expressed as a function of O^, <l>i, '^i, ol, y (Art. 10 (ii)). 

To find L : — Suppose a small arbitrary rotation given to 
the planet about the axis of x, that is, about the line of 
nodes of the plane of maximum areas : then we may repre- 
sent the change in position of the planet by supposing 7 alone 
to vary, 0^, ^^, -^j, a remaining constantf. Thus, the ten- 
dency of L being to diminish 7, we have 

dV 



X = - 



dy' 



the differential coefficient being partial. Hence by the first 
equation of motion, 

da^ 1 dV 
dt • k sin y dy^ 

* The quantity here denoted by h corresponds with what in the Planetary 
Theory is denoted by h, 

f This is equivalent to supposing the plane of maximum areas fixed in the* 
body during the rotation, a. supposition perfectly allowable, since the I'otatioD' 
given to the planet is hypothetical, and has nothing to do with its actual- 
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To find M\ — Suppose a small arbitrary rotation given 
to the planet through an angle h^ about the axis of y, that 
is, about a line through L the centre of gravity of the planet, 
in the plane of maximum areas perpendicular to LM its line 
of nodes : let the eflfect of the rotation be to change the 
position of this plane from MI to mi\ draw Mn perpen- 
dicular to mi. 




Then Mn = 8%, mn = — S'^j, Mm = Sa, the angle Mmn = 7 : 

and we have Z^ = ^^ ^^^ 7* ^'^i = "■ ^^^ cos 7 ; 

rx_ n da 1 dyjr. da, , 

therefore -7- = -. — , -j-^ = — cos 7 j- = — cot 7. 

dx sm7' dx «% 



And 



dx da dx d'^^ dx * 



since we may suppose a and -^^ alone to vary with x > there- 
ore 

,, 1 eZr , dV 

M=—. -, C0t7 j-r- . 

sm7 dx CLYi 

Now (Art. 10 (iii)) t, + (7 =/(A, k, s), 
and 5j, (f)]^ do not involve g ; therefore 

dV dV 

d^i dg ' 

motion. Any other supposition might be adopted with regard to the motion of 
this plane, but this is the most convenient, and will be retained throughout 
this Article. 
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hence -3f = —. j- + cot 7 -j- ; 

sin 7 aa (^ff 

and the second equation of motion gives 

dy _ 1 rfF CO87 dV 
5^ ~ jfclm^ da ■*■ TTim^ ^ • 

To find N: — Suppose a small arbitrary rotation given to 
the planet about the axis of z, that is, about a normal through 
the centre of gravity of the planet to the plane of maximum 
areas : then we may suppose 'yjr^ alone to vary, and since the 
tendency of ^is to diminish 'yjr^, we shall have 

dyfr^ dg ' 
as in the preceding Article ; then the third equation gives 

dk^dV 

dt dg ' 

16. The equation for calculating 7 may also be simply 
obtained as follows : if we refer the motion to the fixed plane 
of reference, and suppose V expressed as a function of 0, if>, 
'^, since k cos 7 is the area conserved on this plane, we shall 
have 

^^(Acos7)=-^. 

Now (Art. 10 (ii)), -^ - a =f{0^, "^j, 7), and 0, ^ do not 
involve a; therefore 

dV^dV 
dsfr da, ' 

the latter diflferential coefl&cient supposing V expressed as a 
function either of 0^, <f>^, -^j, a, 7, or simply of t and the 
elements. Hence 

-ikcosy) = -^^; 
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, « - , dy dV dh 

therefore * ®^^ '^ ^ ~ ^ "*" ^^® '^ "^ 

rfF. dV 



therefore 



dy^ 1 dV cos 7 dV 
di "" A? sin 7 flf a A sin 7 dg' 



The method of this Article is identical with the second 

method given in Art 34 of the Planetary Theory for the de- 

di 
termination of -^ ; we reserve the comparison of the results 

to a subsequent Article (Art. 22). 

17. The following proposition will be useful in finding -^ . 

To shew that jj: ="• 2 -^ , and ^ = — 2 -^ , the ^variables 
being connected together by the equations of Art. 10, (iii). , 



Let u^V^Bh^-^—-^ 

where s = Owg. Then from Arts. 3 and 6, we have 



(1). 



'+'=^(^^)/v|) (2)' 



A-i 
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From equation (2), by dififerentiating under the integral 
sign, 

^ du 
dk 2^^ 'J (ut;)* 

= >J{AB) \ — - — x^ds, by equations (1). 

From equation (3), in like manner, we have 

1/, 8*\{ dv . du\ 
^=-V(^5)J^^-^. j-r^ ds 

, 2uv-(h-pi{Au-Bv) 

-y^^^lik*' -T^ "" 

To simplify this expression, we find from equations (1)^ 

, 8* _ u + v 

jf ^^Au + JBv 
^ "* " A-£ • 

Substituting these values, and reducing, 

1 dl . - 

= 2 ^> from above; 

therefore jt — 2 -^ • 

ak ah 
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Again, from equation (2), we have 

dl 



ds 



V \ wv / 



du 
Diflferentiating equation (4), and substituting for ^ and 

-^ from equations (1), 

also from equation (4), 

'Au\ Au + Bv 



tani^,sec»i^, = y(-g^). 



hence, by substitution, 

dh" 2s/\uvJ' 

therefore 3- = — 2 -^ . 

cLs aft 

18. The following relations between partial differential 
coefficients of the function V, as well as the results of the 

dl 

preceding Ai-ticle, will be required in finding -^ . 

In Art. 10 we have shewn that, in order to express V as 
a function of t and the elements, we may first express it as a 
function of 0^, ^j, -^j, a, 7, by case (ii), and then eliminate 
^i> 0i» '^1 ^y means of the equations in (iii). The first three 
of these equations are 
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If, then, we substitute for 0^^^^^ yjr^ from these equations, 
V will become a function of h, k, a, g^ a, 7 ; and we shall 
have 

/dr\dVd0^ dV d^ dV d^ 
\dhj dff^ dh d<f>^ dh dyfr^ dh 

^dV d(l> dV d^ 
d(f>^ dh dyfr^ dh ' 

since 0^ does not involve h. We have used the brackets to 
distinguish the partial diflferential coefl&cient with respect to 
A, of V expressed as here supposed, from its partial diflfe- 
rential coefl&cient when expressed as a function of t and the 
elements. In order to pass from the expression of K as a 
function of h, k, a, g, a, 7 to its expression as a function of 
t and the elements, we must solve the fourth equation of 
Art. 10 (iii), with respect to a ; thus we shall have 

a=f{h,k, t + l), 

and then 

dV_fdr\ dVda^ 
dh ''[dh)'^ da dh' 

or, writing for ( ^J the value obtained above, 

dh~ da dh'^dif>, dh'^dyfr dh * 
a relation which will be required in the following Article. 

Again, diflferentiating the equation 

a=^f{Kk,t + T) 
on the hypothesis that h and I alone vary, we have 

^'^dh'^dldh* 
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where -^ is the same as would be obtained by direct diflfer* 
entiation of the last equation of Art. 10 (iii). Hence 

da dh da dl dh 

^_dVdl 

dl dh * 

since I is introduced into V only through 8. This relation 
also will be required in the following Article. 

/// 
19. We are now in a position to find-x, and shall for 

this purpose refer the motion to the principal plane of ay. 
We must first, however, obtain an expression for the sum of 
the moments of the disturbing forces about the axis of z. 
Consider, then, F as a function of 0^, <f>^f'^^, a, % and let a 
small arbitrary displacement be given to the planet about 
the axis of z. We may express this by supposing 0^,y^^,a,y 
to remain constant while <f>^ alone varies. Since, then, 
with the usual convention with respect to signs, the ten- 
dency of the moment we are considering is to increase <f)^, it 

dV 
will be expressed by -j-r- . Thus, from Euler's equations of 

motion, 

whence, writing as before, s for Cco^, 

From the results obtained in the case of undisturbed 
motion, we have (Art. 10 (iii)) 
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If we diflferentiate this equation, first considering the 
elements variable, and then considering them invaaiable, and 

equate the results, substituting the above value of -7-, all 

terms not depending upon the disturbing force will dis- 
appear: thus 

dt "" dh dt dk dt ds dff>^ ' 

dh dk 

or, substituting the values of ^ and -^ from Arts. 11 

and 15^ 

^^2— — + —— + — -• 
dt dl dh dg dk dff>^ ds' 



but by Art. 17, 



dk'' dh' da dh' 



also, as in Art. 15, 

dV^^dV 
dg dyfr^' 

and by Art. 18, 

dVd^^^dVd^^ 
dl dh ds dh' 

therefore ^--2('^^* 4. ^^ dyfr dVd<t>A 
therefore^^--2(^^^ + 5;^^^ + ^^J 

ftom which I may be calculated. 

20. To obtain the formula for ^, we might of course 
proceed in the same manner, diflferentiating the third equa- 
tion of Art. 10 (iii) ; but the process would be complicated 
by the fact that '^^ is measured from a moving point and on 
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a moving plane^ and consequently that the form of -^* 

is not the same in the disturbed as in the undisturbed 
motion. We shall therefore employ a different method : 
availing ourselves of the formulae already obtained for calcu- 
lating five out of the six elements, we shall obtain the sixth 
by direct substitution, in the same manner as that in which 
the epoch was obtained in Art. 36 of the Planetary Theory. 

By Art. 10 (iii), we may write 

V^f{t-\'lyg, h,k, a, 7). 
Differentiating, the elements being considered variable, 

d« "" d (« + Z) \ "^ dtj dg dt dh dt'^ dk dt 

dVdoL dVdy 

JL -1- — J 

da dt dy dt * 
differentiating as if the elements were invariable, 

dV^ dV 

dt d{t + l)" 

dV 
Equating the two values of -77 > we have 

dVdl dVdff dVdh dVdk 
dl dt dg dt dh dt dk dt 

^dVd^^dV&i^^ 
dd dt dr^ dt 

dl dh 

Substituting in this equation the values of -^ , -^ , &c., 

obtained in the preceding Articles, we find after reduction 

dg ^ dV cos 7 dV 
dt dk k sin 7 cJy ' 

from which g may be calculated. 
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21. We will here recapitulate the formulaB which have 
been obtained for calculating the elements of the motion : 

,.... dk dV 

r \ ^-. ^^ C0S7 dV 
^ ^ dt dk A; sin 7 d^ ' 

, . da^ 1_ dV 

^ ^ dt ksiny d^ ' 

t '\ <^7 _ ^ ^y cos 7 dV 
^ ' dt A;sin7 da A;sin7 dg' 

22. These six equations agree with those obtained by 
Pontdcoulant in his TMorie AnalyUque du Syst^me du Monde 
(Tome n. p. 189). In Art. 12 we have shewn that the first 
of them can be.expressed in a form in which it is seen to 
be identical with the formula from which the mean distance 
is calculated in the Planetary Theory : it will be instructive 
to compare the remaining equations with those of the 
motion of translation. In order to do this, we shall first 
express the latter in terms of elements having a significa- 
tion analogous to that of the above. Now in the motion of 
rotation h is the element of vis viva, and this in the motion 

of translation is equal to — - : A? is the area conserved on 

a 

the plane of maximum areas, which plane in the motion of 

translation becomes that of the orbit, and the area conserved 

upon it is equal to V/ta (1 — 6)': I is the element added to t 
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in the integration, and this in the motion of translation is 

- (€ — -Bj) : ^ is the element added to -^j, and this corresponds 

to «r in the motion of translation, only we must remember 
that g is measured wholly on the plane of maximum areas 
from its node, while «r is measured on the fixed plane of 
reference as far as the node, and thence on the plane of the 
orbit ; so that «r ~ fl will correspond to g : lastly, a is the 
longitude of the node of the plane of maximum areas, and 
7 its inclination to the plane of reference; in the motion 
of translation the same quantities defining the position of the 
plane of the orbit are denoted by fl and i. In order, then, 
to compare the rotation formulae with those of translation, 
we proceed to replace the elements a, 6, e, «r of the Planetary 
Theory by four new elements A, A;, ?, g^ connected with the 
former by the relations 

a 
A* = /ia(l-e'), 

Z = - (€ — -Bj), 

Let jR' express the form which the function R takes 
when the new elements are substituted for the old: then 
we have 

dR dR dh , dE dk . dR dl ^ 
da dh da dk da dl da 



~'i?~dh^'^ dk"^ 2na dl * 



an 



* In -T- , which occurs only in the formula for the epoch, the differen- 



tial coefficient is supposed to be taken with respect to a only in so far as o 
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dR 


dR' 

dk 


dk 
de~ 


na*e 


dR' 


de 


Jl-e* 


dk* 


dR 


dR' 


dl 


1 dR' 




de ' 


^ dl 


de~ 


n dl ' 




dR 

dvr 


dR! 

~ dl 


dl 

der 


dB dg 
dg dvr 




1 


IdR' 
~ n dl 


dR: 
dg* 




dR 

da' 


dR' 

'da 


dB! 


dg 

da 






dR 


ds: 







dH dg * 

Now if we diflferentiate the above expressions for A, k,l, g, 

and then substitute the values of -^ , -rr, &c. obtained in 

dt dt 

the Planetary Theory, and those obtained here for -r- , -j- 5 
&c., they will reduce to 

dJt A^^ 

di^ 'dl' 

dl ^ o ^^' 
di'"^ cLh' 

dk^dR 
dt dg ' 

dg dR' cost dR' 
dt" dh k sin i di ' 

occurs explicitly mR: if we suppose a to vary also as contained implicitly in n, 
this difEerential coefficient will include the term proportional to the time, which 
may therefore with this understanding be omitted. In Art. 37 of the Platietary 
Theory this term was removed by a different transformation. 



30 THE eakth's motion of rotation. 

dCl^ 1 dR 
dt k sin i di ' 

di ^ 1 dR' cost dR 
dt~ h sin { dOi h sin i dg 

On comparing these with the rotation formulae, it will 
be seen that they are identical, with the exception only that 
the sign of a diflfers from that of n ; and this is accounted 
for by the fact that a and n are measured in opposite direc- 
tions. Thus, by employing elements having a like significa- 
tion in the two motions of translation and rotation, we have 
arrived at the very remarkable result that the complete solu- 
tion of the problem of Planetary perturbation, whether in the 
motion of translation or of rotation, is expressed by the above 
simple formuloe. 



PARI IL 

APPLICATION OF PRECEDING RESULTS. 
PRECESSION AND NUTATION. 

23. The formulae obtained in the first part are suf- 
ficient completely to determine the motion of a planet or 
other rigid body about its centre of gravity. They are per- 
fectly rigorous, subject only to the hypothesis that the dis- 
turbing bodies may be supposed to attract as if condensed into 
their respective centres of gravity ; a hypothesis admissible if 
these bodies are either very distant or nearly spherical in 
form. But though the formulae are exact, they can be in- 
tegrated only by approximation. We propose, therefore, in 
the second part to restrict ourselves to the particular case of 
the Earth, taking advantage of such of the results of obser- 
vation as may be required to enable us to approximate. In 
order to treat the problem fully we shall consider, first, the 
motion of the axis of rotation in the Earth itself, with the 
velocity of rotation about it ; secondly, the motion of this 
axis in space. The first is of special interest, since any 
change in the position of the axis in the Earth, were such 
change possible, would aflfect the permanence of terrestrial 
latitudes ; any change in the velocity of rotation would aflfect 
the length of the day. The second is of great importance to 
astronomers, since it establishes the fact that the first point 
of Aries, or vernal equinox, to which they are accustomed to 
refer celestial longitudes, is not a fixed point. 
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Stability of ike axis of rotation in the Earth and 

of the velocity about it. 

24. We proceed then, first, to consider the motion of 
the axis of rotation within the Earth, and shall be able to 
shew that, in so far as it depends upon the attractions of 
other bodies, this axis can never separate appreciably from 
the axis of figure, and that the velocity about it must always 
remain appreciably constant. 

We have from Ai-t. 3 

Eliminating to^ from these equations, 

Now if (7 be either the greatest or least principal mo^ 
ment, both the terms of the left-hand member of this 
equation must always retain the same sign ; and if we take 
C about the Earth's axis of figure it will be the greatest and 
this sign will be positive. We may therefore write 

thus ^(0--4)o>,* + 5(a-5)ft>/=((7--4)((7-5)e'. 
It follows from this that 



(Wj can never be greater than e ,k/ : 



G'-B 



coj can never be greater than e 






If we neglect the* disturbing force, A, h and therefore e 
are constant, and the value of e will be found by substituting. 
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the values of a>j, o>, as determined by obsenration at any 
given epoch. Now the most delicate observations have 
hitherto shewn no appreciable separation of the axis of 
i*6tiation from the axis of figure : hetice Op o), are at present 
inappreciable; thus e is so, and we conclude that, indepen- 
dently of the disturbing force, w^, o), must always remain 
inappreciable. 

If in the equations of Art. 6, we write to^ = 0, o), = 0, we 
have also ^^ = ; so that, if the disturbing force be neglected, 
we may consider the plane of maximum areas as coincident 
with the Earth's equator. 

c 

25. Let us now examine the effect of the disturbing 
force upon the angular velocities w^, ©,. By the preceding 
Article 

((7--4)(C-5)£»=CA-*». 

Differentiating, and substituting the expressions for -^- 

dk 
and -TT obtained in ihe first part (Art. 21), 

If the approximation be carried to the first power of the 
disturbing force, we may, in calculating the second member 
of this equation, suppose ^^ = : thus the equations of Art. 7 
give 

Also, if tt be the mean angular velocity, we may to the 
same order of ajqproximation suppose m^^ni thus, from the 
equations (A) of Art. 5, ^ and are constant and 

C. 3 



84 THE earth's motion op rotation, 

whence by integration 

e being the value of ^ at the epoch from which t is measured 
Again, by Art. 10 (iii), 

t + I =f{h, k, s). 

dV^ dV ^dV^l d{V) • 
dff dy^^ d4> n dt * 

Now (Art. 3) A» = ^ V + -» V + C^ V 

to the same order of approximation : therefore 

k=Gn. 

„ ^dV .dV ^d(V) kd{V) ^ 

•Hence C -jf — k-j- — C ,, ^ V-^ = ; 

dl dg dt n at 

and therefore ^r = : 

dt 

whence it follows that to the first order of the disturbing force 
e continues absolutely constant. 

When the approximation is carried to the order of the 
square of the disturbing force*, it is found that e* can con- 
tain no term proportional to the time, nor any inequality, 
either secular or periodic, raised to the first order in the 
process of integration. Hence we may consider 6* as prac- 
tically constant, and thus conclude, «^s in the preceding 
Article, that o^, ©^ will always remain inappreciable. 

* See PentdcoQlant*B Syttime du Monde, tome ii. p. iiS. 

/J 
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It follows, as in the preceding Article, that to the same 
degree of 9.pproximation the plane of maximum areas may 
be regarded as coincident with the Earth's equator. 

26. We have now to examine the effect of the disturb- 
ing force upon the velocity of rotation. Since g)^, g), are in- 
appreciable, this velocity will be equal to ©g ; thus, as in thq 
preceding Article, we have 

k 

therefore ^» = ^ - = ^^ ^ (Art. 21) 

^ 1 diV) 
Cn dt • 

Now considering F as a function of 0, <j>, '^, if our ap- 
proximation be carried to the first power of the disturbing 
force, we have 

= y, ^ = nt + c, '^ = a; 

so that V will involve t only under the form nt + c. It fol- 
lows that — TT- will consist of a series of terms of the form 
dt 

sm ^ ^^ 

where p is some positive integer. And this term will give 

27r 
rise to inequalities whose period is — , that is, a day or 

fi'action of a day. Now the most delicate observations have 
hitherto failed to detect any such inequality : hence we 
must altogether reject such terms, and conclude that G)g will 
always remain sensibly constant. 

27. The preceding Articles prove the stability of the 
axis of rotation in the Earth and of the velocity about it, so 
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far as the motion depends upon the attractions of distant 
bodies ; but there is another disturbing influence to be found 
in the friotion produced by the tides* The effect of this (see 
Thomson's and Tait's Natural Philosophy) would be to retard 
the Earth's velocity ; but whether to an appreciable extent 
is yet uncertain. The fact, discovered by Professor Adams, 
that a portion of the acceleration of the Moon's mean motion 
is yet unaccounted for, led Delaunay to suggest that this poi"* 
tion might be apparent only, and really due to the retarding 
effect of tidal friction upon the Earth's velocity. But other 
possible explanations have been given, and to pursue the 
subject further would be beyond the scope of the present 
treatise. 



Precession and Nutation, 

28. We have now to consider the motion of the Earth's 
axis in space, resulting in the phenomena of Precession and 
Nutation. We have seen that, to the first order of the dis- 
turbing force, the plane of maximum areas may be con- 
sidered as coincident with the Earth's equator, and that 
consequently to this order of approximation we may write 

tf = 7, -^ = a. 

Hence the motion of the equator may be calculated by 

formulae (v) and (vi) of Art. 21. Now it may be shewn as in 

Art. 25 that 

dV_ld(V) 

dg n dt ^ 

ftud may therefore be neglected, since (as has been proved in 
Ai-t. S6) it can give rise only to terms whose period is a day 
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or fraction of a day, which are wholly iusepsible. If, then, 
we write y^ and fpr a and 7 respectively, these formulae 
become 

df 1 dV 

dt ^ Asind d0* 

■ ^^ 1 dV ^ 
dt ibsin^ rf^' 

or, since to the same order of approximation h = CSw, = CW, 

d'^^ 1 AV 

dt "" Cn^me de ' 

d0^ 1 dV 

dt" On sin ^ dy^ ' 

39. Thes^ elegant equations, due to Poissou^ giye '^ 
and 0y and thus determine the motion of the Earth's axis in 
space. Before we can integrate them it Mrill be necessary to 
obtain an expression for V. 

Let then m be the ipass of the Earth, ni that of the dis*- 
turbing body which we shall suppose to be the Sun ; r the 
distance between the centres of gravity of these bodies, r the 
distance from the Earth's centre of gravity of an element Sm 
of its mass, p the distance of the same element from the 
centre of gravity of the Sun : also let the inclination of r to 
r' be denoted by v. Then by Art, 9, 

P 



(r* — 2r r cos V -f r^* 
= ^5|8/n(x-.yCosr + ^?) |. 
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Now since / is very large in comparison of r, we shall 
neglect all powers of - above the second : thus 

F=5 SJSm (l + ^ cost;- I ^3 + ? ^cos»t;)| 



= — S(S»»)+-;^2(Sm»rcosv) -g-TaS (8m. r') 



+ 2^ 2 (Sw . r^ cos'v) 

ft f 

= — 2 (Sm) + -75 2 (S»» . r cos i?) + —, 2 (Sm . r') 

— rt-75 2 (Sm . r' sin*i;). 

Of these terms -T2(Sm) and -^2(Sm.r') may be omitted, 

since they cannot contain the angled 6^ <f> or yfr, and would 
not be affected by any arbitrary rotation such as is supposed 
in Art. 9 : also since r is measured from the Earth's centre of 
gravity, 2 (Sm . r cos v) = 0. Hence we may write 

F=-|^2(Sm.7^sin*t;) 

if Q denote the moment of inertia of the Earth about the 
line joining its centre of gravity with that of the Sun. We 
see, then, that the effective portion of the disturbing func- 
tion is proportional to this moment of inertia. 

30. We have hitherto made no assumption respecting 
the Earth's figure. We learn from pendulum observations 
and geodetic measurements that it is approximately an 
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oblate spheroid; the moments of inertia about all axes in the 
plane of the equator are therefore nearly equal, and we may 
take^^«B^. If^ then, £ be the Sun's declination, 

C = .4cos'S + (7sin'S: 

also if n' denote the Sun s mean motion, we have approxi- 
mately 

therefore — rs — = W\ 

or, neglecting the Earth's mass in comparison with that of 
the Sun, 

Hence, by substitution, 

F=,«^ (^ cos'S + (7sin«S). 

31. By Art. 9) the moment of the Sun's disturbing force 
about an equatorial diameter of the Earth perpendicular to 
the line joining the centres of gravity of the two bodies, and 
tending to increase S 

It is clear that the moment about any axis perpendicular 
to this is zero, since an hypothetical rotation about any such 
axis could produce no change in S. If, then, the Sun's dis- 
turbing forcer be Teduced to a parallel force through the 
Earth's centre of gravity, and a couple, the tendency of the 
latter will be to give to the equator a rotation towards the 
ecliptic about an equatorial diameter perpendiciular to the 
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line joining the centres of gravity of the two bodies. The 
moment of this disturbing couple is 

3n» 



2 



{C-A)sm2S. 



32. We shall now return to the equations of Art. 28, 
and deduce from them the motion of the Earth's axis in 
space. The results obtained will be of sufficient accuracy for 
our present purpose if we suppose the ecliptic, to which plane 
the motion will be referred, to remain fixed. 

Let yfr denote the longitude of the first point of Aries 
measured from some fixed point on the ecliptic, X the Sun's 
longitude measured from the same point, the obliquity of 
the ecliptic. Then from the S{dierical triangle formed by the 
intersection of the equator, the ecliptic, and a declination 
circle through the Sun, we have 

sin S = sin ('^ "^ X) sin ^ ; 
therefore cos S -1^^ = sin (i|r — X) cos ^, 

cos5 3-r = cos (-^-^Xjsii^ft 

If -^ — X = Z, then I will be the Sun's longitude measured 
from the first point of Aries in the opposite direction to that 
in which yfr and X are measured : thus 

cos 5 ^ as sin ? cos ^, cos 8 ^r =» cos Z sin 3. 
Hence -3z^^ ^=s — 8n''((7— -4)sin8sin?cos5, 
-5-r- = -js- -TT =-"3n"(<7— -4)sin8cosZsin^» 
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Substituting these values in the equations of Art. 28, 

dytr 8n'" C-A . ^ . , ^/, 
-T" =?= — . — 7y — sm 5 sin f cot a 
at n V 

Sn C — ^ . - - ^ 
« — . — p.— sm' e cos r 
n C 

8«'' C—A y, ,- ^J. 

~ "o— • — 7r~ cos ^ (1 — cos 2i}. 
de 3a' O-^A . . , 

J- = . 7^ — smOCOS t 

5n? C-A . , , . ^ 

= . — 71 — sm 6 cos I sm ^ 

w G 

= — -7i— . — TT- sm ^ sm 2fc. 

Let / denote the mean value of the obliquity, then since 

n C — A 

both — and — —^ — are very small quantities, we may in the 
n G 

second members of these equations write 

thus -^ = ^— . — Yr~ cos 1(1 — cos zn t), 

dd 8n'» C-A . r ' o u 

By integration, 

^ ^ Sn' G^A . r 
= 1 + -2 — . — PS — . sm I, cos zn t 
4sn U 

y being measured from the position of the first point of Aries 
at the epoch from which the time is reckoned. 
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33. These formulae determine the motion of the Earth's 
axis in space. Although the method by which they have 
been obtained has the advantage of connecting the two prob- 
lems of translation and rotation, it may be worth while to 
arrive at them by an independent process. 

Since the Earth is nearly an oblate spheroid, we shall, as 
before, take C for the greatest principal moment, and sujp- 
pose 5 = -4, so that all axes in the plane of the equator are 
principal axes. If, then, taking the plane of the equator for 
that of ajy, we suppose the axes of x and y to revolve with an 
ai^ular velocity 0^ about the axis of «, we have the equations 
of motion* 



dt 



0^* = N. 
at 



Now if we take for the axis of x the projection of the 
Sun's radius vector on the plane of the equator, we have 

i = 0, j\^=0, 
Jlf=-^(C7--4)sin2S, (Art. 31); 

and since 0^ is very small and occurs only in connexion with 
<»j i^rid a>, which are also very small (Arts. 24 and 25), we 
may write 

^, = n its mean value. 

* See Riouth's Rigid Dynamics, Art. 107 : the equations there given reduce 

dx 
to the above if ^s is written for «&+ -^ . 
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Thus the third equation gives 

whence ©, = const. = n suppose ; 

and the first two become 

6r, denoting ^ by p, 

C—A . 
Now — ^T — • is a very small fraction ; n , the Sun's mean 

motion, is very small in comparison of n ; and 8 vaiies very 
slowly; we shall, therefore, in integrating these equations, 
take no account of its variation. Thus, eliminating g)^, 

the integral of which is 

a, = A cos (pt — B), 

where A and B are constants of integration. This term is 
independent of the disturbing force; but having for its period 

— , which does not differ much from a day, it must be re- 
P ■ 
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jected, since, as we have already remarked in Art. 26, such 
terms are altogether insensible. Hence we conclude that 

and by substitution in the second equation. 



©J = -Tz — . — : — sin 2S 



2p' 


■ A 


Sn" 


C-A 


2» • 


C 


^H* 


a- A 






sin2S 



-2n- C ^^^' 

powers of n' above the second being neglected. We see^ 
then, that the effect of the Sun's attraction is to cause the 
Earth's axis to travel in a plane perpendicular to that of the 
disturbing couple which it produces, and with an angular 
velocity proportional to its moment. This velocity, though 
so Bmall as to be insensible to the most delicate observations, 
yet leads to values of ^ and 9 which can by no means be 
neglected, 

34. We proceed to determine the motion of the axes in 
apace. Substituting in the equations of Art. 5, we bavQ 

-±-Binu^ — • — 71— . sin o cosdsm6, 
at ^ n ^ 

dd 3n' G-A . ^ . . 

►1- = . ^ ^u . sm cos o cos o. 

at n V ^ 

Now from the spherical triangle formed by the interseC'^ 
tion of the equator, the ecliptic and a declination circle 
through the Sun, it is easily seen that 

sin ^ == tan h cot 6^ 
sin S = sin Z sin 0, 
cos Z = cos ^ cos S, 
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By means of these equations we obtain 



dt n 



. cos sin' Z, 



whence, as in Art. 32, 

, 3»'» C-A J. , 3ft' G-A r • o '. 

^ = /+-J — • — -; — sin /cos 2/1 1 
4n 

35. It appears from these formulaB that the motion of 
the Earth's axis is of two kinds ; partly secular, partly 
periodic. It is convenient to consider these sepai'ately. The 
former affects the equinoxes alone, and, being proportional to 
the time, indicates uniform motion; this motion is one of 
regression, since -^ has been measured in a direction contrary 
to that of the apparent motion of the Sun. Considered with 
reference to the apparent diurnal motion of the stars, the 
effect is to place the equinoxes in advance of the position 
they would occupy if fixed : hence it obtained the name of 
the Solar Precession of the Equinoxes. The latter furnishes 
corrections both on y^ and 0, which go through all their 
changes in half a year : it is called the Solar Nutation of the 
Earth's axis; the correction on '^ forming the Nutation in 
Longitude, that on the Nutation in Latitude, 

36. We will now examine the effect produced by the 
Moon's action on the motion of the Earth*s axis. Since the 
investigations which have been given of the effect of the 
Sun's disturbing force contain nothing to restrict their gene- 
rality except the special assumptions made with regard to 
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small quantitiesi we shall first consider what modifications 
are required in order that they may be applied in the case of 
the Moon, 

Let '^\ ffy ri\ T denote relatively to the plane of the 
Moon's orbit the same quantities which relatively to the 
ecliptic have been denoted by -^^ 0, n\ I\ also let m, rri' be 
the masses of the Earth and Moon : then^ as in the case of 
the Sun, we have approximately, 

2ir 27rr"^ 



but we cannot in this case neglect tti, as it is in fact much 
larger than w" : retaining it, we have 






if m=^\m". Also n" is small, though not so small as w' in 
comparison of n. Hence to detennine the motion of the 
Earth's axis with reference to the plane of the Moon's orbit, 
we have 

Xt = TT — 7= r-. , 7y — . COS 1 t 

^ 2n (1 + \) C 

on C — jA. j.f ' c% 't, 

-j — . i sin ZTi t, 



4n(l + X.)' 

^ = 1 +ir-7^ — ^x- — 77— •smi cos2n t. 
4n (1 + \) (J 

The periodical terms in these equations go through all 
their values in half a month, and are so small that they 
are usually neglected. Thus we may consider the incli- 
nation of the plane of the Earth's orbit to that of the Moon 
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jas constant and equal to its mean value, and the precession 
as uniform and given by 

whence, the velocity of precession 

dt ~"2n(l + \)- G •'^^^^' 

These results, deduced from the corresponding formulae 
in the case of the Sun, suppose the plane of the Moon's orbit 
fixed Now the line of nodes moves too rapidly to allow of 
this hypothesis, but the only efiect of considering its motion 

would be to add to --j- a term depending upon its velocity 

and not upon the disturbing force of the Moon upon the 
Earth. Since our object is to trace the eflFects of this force 
only, such terms must be omitted. 

37. It now remains to determine the motion of the 
Earth's axis with reference to the ecliptic. 

Let be the obliquity of the ecliptic, -^ the longitude 
of the equinox, a the longitude of the node of the Moon's 
orbit measured from the same origin as '^^ i the inclination 
of the Moon's orbit to the ecliptic, /' its inclination to the 
equator. Then supposing '^' measured from the node of the 
Moon's orbit on the ecliptic, we have by Spherical Trigono- 
metry (as in Art. 7) 

sin (i|r — a) sin ^ = sin J' sin i|r' (1), 

cos = cos i cos/' ■- sin i sin/' cos '^' (2), 

cos /' = cos cos / + sin sin I'cos (-^ — a) (3). 

"We must now difierentiate these equations in order to 
express -^ and -^ in terms of -^ . In so doing we. shall 
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omit the terms involving -r and -7-, since these quantities 

depend not upon the disturbing force, but upon the motion 
of the Moon's orbit. Thus 

cos (-^ — a) sin 6 -^ + sin {^ — a) cos ^ -7- 

sssini'cos'^' -^, 

sm ^ ^- = — sm t sm I sm ylr —j- . 
at ^ at 

Writing XI for -^ — a, we obtain 
cos 12 sin ^ -J- = (sin /' cos yjr 



+ sin i sin F sin '^' cot sin XI) -^ 

tZv 

= (sin /' cos yfr + sin z cos sin' X2) -^ by (1), 

'cos i cos /' — cos 



sinz 



+ sinz cos5 — sin/ costf cos'Xlj -^ by (2), 

/cos t cos T' — cos ^ cos'/ . . ^ , ^N e?>t' 

= ( ; — ; sm t cos ^ cos Xi —j- 

\ ' smt / dt 

= (cos/sintf cosXl — sin/costf cos'Xl)-^ , by (3) ; 

therefore -^ = (cos / — sin / cot cos XI) -^ . 

., d0 sin / sin J' sin yj/ d^jr' 

dt ^ sin dt 

= — sin/sin fl -^ by (1). 

1/ f ' 
We must now substitute the value of -y- : thus 
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-T- =s— t: — ^T- — TT- .cosi^ (cost — Sin t cot^cosllK 
dt 2n(H-X) ^ ' 

It remains to eliminate jT. Now 
cos T (cos t — sin * cot S cos H) 
= (cos 5 cos t + sin sin t cos ft) (cos t — sin z cot cos fl) 

= cos cos** — s -^ — ?r sin 2/ cos 11 — cos ^ sin*t cos*fl 

2 smd 

af !• 1 • « A 1 cos2^ . gj, ^ 
s= cos ^ cos'» — H sm" * } — 5 — ; — ^ sm 2* cos H 
\ 2 / 2 sm& 

— ^ cos 5 sin"* cos 2fl. 

Also cos/'sint sinfl 
= (cos cos t + sin sinz cos A) sin i sinO 

1 ... 1 . . . . 

= ^ cos sin 2i sin fi + ^ sin sin'* sin 211. 

^"^^ ■^=2jr(rTx)--o--^"^^(^«* 

1 . . 1 

-^ 5 sin"* — cot 20 sin 2* cos fl — ^ sin*» cos 2fl), 

-^ = — "1 — 71 — TT • — 77 — • (cos sm 2* sm il 
at 4w(l + \) C/ ^ 

+ sin d sin'* sin 211). 

Let V be the mean retrograde velocity of the Moon's 
nodes, H^, the value of fl at the epoch, then we may write 
il^vt + £l^; also for we may write its mean value I; and 
since * is a veiy small angle we may neglect sin'* when mul- 
tiplied by a periodical term. Thus 
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dir 3n"" C^A J.. ,. 1 . ,. 

-T- = -^ — 7^ — ^ • — 7i — • COS 1 {COS t — 5 sm'e 
dt 2n{^-^\) C ^ 2 

— cot 27 sin 2^ cos (j/< + flj}, 
d0 3> C-A r • o- • /. ON 

The integrals of these equations are 

Sn"* C-A j.f, ,. 1 . ,^, 
^ = 2>.(1 + X) -~C~ •<'os/{(cost-gsm»0< 

cot 27 sin 2^ sin {vt + 11)}+ const., 

= I+-. — 7= — r^ . —7y — cos7sin 2^ cos (vt + flj. 
4n(lH-\) Cv ^ ®' 

38. These equsftions determine the motion of the Earth's 
axis due to the attraction of the Moon. They are similar 
to the equations expressing the Sun's action, and the 
remarks made in Art. 35 might, mutatis mutandis^ be re- 
peated here. The Lunar Precession is 

3n'" G-A j( ,. 1 . ,A, 
-ET-fk — ^^« — ?? — . cos / I cos'* — 5 sm'i t. 
2n(l + X) (7 V 2 / 

K we add to this the Solar Precession (Art. 32), we find 
for the whole permanent eflfect of the Sun and Moon upon 
the equinoxes 

3n G^A ^( n"' ( ,. 1 . , .\ . n\^ 

This is called the Luni-solar Precession, to distinguish it 
from the Precession due to the secular motion of the ecliptic 
in consequence of the attractions of the Planets, and which 
is called Planetary Precession, 

Adding together the eflfects of Lunar and Solar Nutation, 
we find for the whole Nutation in Longitude 
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3 C-A .f n"' 



. cos/' 



2* C {np{l + \) 



cot 27sin2i sinfl 



where II and denote respectively the mean longitude of the 
ascending node of the Moon's orbit and of the Sun. 

Similarly, the whole Nutati(jri in Latitude is 
3 C^A ^( n'" . ^. ^ n' 



4* G 



• cos I \ — 7^ — r-r sin 2^ COS ft H coS 2 >• . 

[nv (1 + A-) n ) 



In each of these expressions the first term, due to the 
action of the Moon, is the most important, since n^' is larger 
than both n' and v. 

39. If we consider only the motion of the Earth's axis 
due to Precession, it appears from the preceding formulae 
that it maintains a constant inclination to the pole of the 
ecliptic, and describes a right circular cone about it with 
uniform velocity. An axis possessing this motion exactly we 
shall term the mean axis of the Earth. 

The motion of the axis due to Lunar Nutation can now 
be exhibited as follows : — 

Let aj = — y. — 7^ — TT. — 7y — . COS 2/ sin 2t sin ft, 
4 nv{l + X) (J 

V = T . — 7^ — :;-v . — Ti — • COS 1 sm 2i COS H : 
^ 4 wz/ (1 + \) C 

then, taking unity as the radius of the celestial sphere, x 
and y are the small linear spaces traversed by the intersec- 
tion of the Elarth's axis with the circumference in two direc- 
tions at right angles. Eliminating II from these equations, 
we have 

a? , y" ^ fS n!'' C^A . .V 
c"^?2l"^^5?7""Uw(l-i-\)' -^'^^V 



52 THE earth's motion OF ROTATION, 

Hence, in consequence of Lunar Nutation, the extremity 
of the axis may be considered to move in an ellipse whose 
semi-axes are in the ratio of cos 21 to cos /. The centre of 
this ellipse is at the point of intersection of this mean axis 
with the celestial sphere, and its plane a tangent to the 
sphere at that point. By supposing the mean axis to de- 
scribe the cone uniformly, while the true axis describes this 
ellipse about it, the real motion will be represented. This 
conception is due to Bradley, who arrived at it by obser- 
vation. 

40. The annual value of the Luni-solar precession (see 
Art. 38) 



2n' a 



cos/-^~27= — r-v cos'« — s sm't ^--rt-n. 



ObseiTation gives about 50".l as the numerical value of 

this expression. Hence, by substituting the known values of 

C — A 
n, n\ V, J, i, we have a relation between — ^y — ^^^ ^> ^7 

means of which either may be determined when the other 
is known. 

41. We have taken no account in our calculations of 
the Precession caused by the attractions of the planets on 
the Earth, since it is too trifling to be appreciable. 



THE END. 
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